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Abstract 

I review the most recent progresses in the calculation of four-fermion processes in e + e~ collisions. 



Introduction 

Final state four-fermion processes represent an important ingredient when studying high energy e 



+ P - 



collisions. They enter in the analysis of the Z-peak observables at LEP, as (mainly) QED pair corrections 
to 2-fermion processes. Such a set of contributions, together with the full set of clcctrowcak (EW) loop 
corrections, allows very stringent tests the Standard Model (SM) at the radiative level. 

At LEP2 energies, all relevant signatures such as WW, ZZ, single-VF, single- Z and Higgs production 
manifest themselves as four-fermion final states. These same processes are also relevant for precision 
measurements at the Linear Collider (LC) and as a SM background to searches. 

In the following, based on the CERN Report of ref. [[I], I review the most recent achievements in the 
computation of four-fermion processes. Table [l] contains a list of contributing codes. Emphasis is put 
here on the available tools for estimating the theoretical errors to be associated with the four-fermion 
observables, in view of the final LEP analysis, but also on the improvements needed at the LC. 

The first two sections are devoted to VF-pair and single-VF signatures, while, in the last two parts, I 
cover topics on four-fermion production plus 1 visible photon and Z-pair final states. 



Code 



BBC | 
CompHEP 0] 
GENTLE | ] 
GRACE g 
KDRALW/YFSWW 
NEXTCALIBUR 
PHEGAS/HELAC 
RACDDNWW U 

SWAP |[o) 

WPHACT 
WRAP JlO 
WTO p 
YFSZZ 
ZZTD [|l 
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Table 1: Contributing programs. 



1 VK-pair production 

When collecting data at y/s = 189 GeV, the LEP2 collaborations observed a deficit in the number of 
events, with respect to the SM predictions. This fact triggered a re-analysis of the available tools for 
calculating the total cross section oww ■ At that time, a theoretical 2% error band was assigned to this 
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observable, two times bigger than the experimental error. The estimate of the theoretical error was based 
on the GENTLE/4FAN inclusion of QED Initial State Radiation (ISR), without any attempt to take EW 
contributions into account. Therefore, it was immediately clear that the computation of the genuine EW 
effects was needed to match the experimental accuracy. On the other hand a full four-fermion one-loop 
EW calculation was (and still is) beyond reach, and including only the WW-like diagrams violates gauge 
invariance. The solution to this problem is represented by the so called Double Pole Approximation 
(DPA) The DPA isolates the poles at the complex squared masses, with gauge invariant residues 

which are then projected onto the respective on-shell gauge invariant counterparts. The projection is 
from the off-shell phase space to the on-shell phase space. Even though such a procedure is strictly gauge 
invariant, the projection procedure is not unique. However, the ambiguity is small, namely 

For example, in the case of a single unstable particle, the fully re-summed amplitude can be rewritten 
as follows 
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(1) 



where M and M are the bare mass and the complex pole of the instable particle, and Z the wave-function 
factor. The first term is the gauge invariant single-pole residue (on-shell production and decay of the 
unstable particle). The second term has no pole and can be in principle expanded in powers of p 2 — M . 

Applying DPA to W-pair production means that only the double-pole residues of the two resonances 
are considered, and one-loop EW contributions included there, for which only (available) on-shell cor- 
rections are needed. The corrections to be included fall in two different classes, namely factorizable 
contributions, in which the production, propagation and decay steps are clearly separated, and non- 
factorizable contributions, in which a photon with energy < is emitted (see figure |l|). The DPA is 



(a) 




production 



decays 



(b) 




Figure 1: Factorizable (a) and non factorizable (b) contributions to VF-pair production, 
not reliable at the VF-pair threshold, where the background diagrams get important. The expected DPA 
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uncertainty above threshold is of the order O (f-;p^7 m (' ' ')) < 0.5%, in fact, when > 2Jl% + nTw 



0.1%. 



with n — 0(3 — 5), the background diagrams are of the order ~ 2 m ^ n (' ' 

Very far away from resonance, the DPA cannot be used any more. 

At LEP2 energies, the inclusion of the DPA formalism in RAC00NWW, BBC and YFSWW allows to lower the 
theoretical uncertainty on <tww from 2% to 0.5 %, in much better agreement with the data. In figures ^| 
and [jj and table ^ we show examples of comparisons among the codes. 
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Figure 2: Relative 0(a) corrections to er(e + e~ — * f M /x + T _ P t ). 



no cuts 


c t ot [fb] 


final state 


program 


Born 


best 




YFSWW3 


219.770(23) 


199.995(62) 




RacoonWW 


219.836(40) 


199.551(46) 




(Y-R)/Y 


-0.03(2)% 


0.22(4)% 




YFSWW3 


659.64(07) 


622.71(19) 


udfi~ ' 


RacoonWW 


659.51(12) 


621.06(14) 




(Y-R)/Y 


0.02(2)% 


0.27(4)% 




YFSWW3 


1978.18(21) 


1937.40(61) 


udsc 


RacoonWW 


1978.53(36) 


1932.20(44) 




(Y-R)/Y 


-0.02(2)% 


0.27(4)% 



Table 2: Total cross section at y/s = 200 GeV without cuts. 

In conclusion, with the help of the DPA, a theoretical accuracy at the level of 0.5 % on aww is 
reached, as required by the LEP2 collaborations [l]]. The error decreases with increasing energy, giving 
the following estimates of the theoretical uncertainty on ctww 

0.4 % at y/s = 200 GeV, 0.5 % at y/s = 180 GeV, 0.7 % at y/s = 170 GeV. 

A theoretical uncertainty of the order of 1 % must be assigned to the distributions. 

A full four-fermion one-loop EW calculation is still missing, but it is required for high precision 
measurements at the LC. 

2 Single-W 7 production 

The experimental single- W signature is better defined with the help of table pi The contributing Feynman 
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Figure 3: Distribution in the W + invariant mass at y/s = 200 GeV for the ud[i final state. 



Process 


diagrams 


cuts 


eevv 


i-channel only 


E(e+) > 20 GeV, |cos0(e+)| < 0.95 


evixv 


i-channel only 


E(H+) > 20 GeV 


evTv 


i-channel only 


E ( t+) > 20 GeV 


evud 


i-channel only 


M(ud) > AbGeV 


eves 


i-channel only 


M(cs) > AbGeV 



Table 3: Possible single- W processes. |cos0(e )| > 0.95. 




s-channel t-channel 

Figure 4: s-channel and i-channel diagrams in single-W^ production. 

Only the t-channel contribution is included in the definition of the single-VF processes. This set is 
explicitly drawn, for the evud final state, in figure [|. The reason of this diagram based definition is that 
it allows an easy combination of different processes from different experiments. Notice that, in order to 
preserve gauge invariance, all i-channel diagrams must be included, not only those where a W boson is 
produced. 

Since in the limit of vanishing electron masses the t-channel diagrams blow up in the forward region, 
the computation of single- W processes is a challenge from a technical point of view. In addition, including 
by hand a width for the W boson breaks the U(l) gauge invariance, so that gauge preserving schemes 
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must be applied, such as the Fixed Width (FW) and the Imaginary Fermion Loop (IFL) approaches |l6| ]. 

Generally speaking, the technical questions are well under control, both for total cross sections and 
distributions, as shown by the comparisons in tables and figure |^. 

As for a realistic modelling of this process, the main difficulties are due to the presence of different 
scales. Usually, when studying high energy processes, part of the higher order corrections can be reab- 
sorbed in the Born approximation by using the so-called Gp scheme. In such a scheme Gp, Mz and Mw 
are input parameters, while the weak mixing angle and apM are derived quantities: 

s 2 w = 1 - Ml I Ml , a EM = V2 G ? M w*w . (2) 

In the presence of low i-channel scales such an approach fails, since the choice a(t ~ 0) ~ 1/ 137 is certainly 
more appropriate for the i-channel diagrams of figure |4| The question is therefore how to consistently 
include the running of aEM in single- W four- fermion processes, without breaking gauge invariance. 

An exact and field-theoretically consistent solution is represented by the Exact Fermion-Loop (EFL) 
approach of ref. fll7| , where the whole, gauge invariant set of fermion one- loop corrections is taken into 
account, including the vertices. 

In table || the results of FW and EFL are compared for different cuts on the outgoing electron. 

There are also approximate solutions, such as IFL Q jl8) - where a(t) is used for the i-component 
and OiG F for the s-component - or the Modified Fermion Loop (MFL) approach by NEXTCALIBUR, that 
only includes the leading self-energy like corrections from the EFL plus an effective vertex to preserve 
the f/(l) gauge invariance. 

IFL Q and MFL coincide numerically at LEP2 energies, but slightly disagree with respect to EFL (up 
to 2% at LEP2). Comparisons at different energies are presented in tables [?] and ||. 
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s/s = 183 GeV 


s/s = 189 GeV 


y/a = 200 GeV 


NEXTCALIBUR 


26.483 ±0.041 


29.679 ±0.047 


35.893 ±0.048 


SWAP 


26.47 ±0.04 


29.70 ±0.04 


35.93 ±0.05 



Table 4: Cross-sections [fb] for e + e — » e u e ^ + f^. 





s/s = 183 GeV 


s/s = 189 GeV 


s/s = 200 GeV 


NEXTCALIBUR 


26.422 ±0.035 


29.655 ±0.046 


35.954 ±0.052 


SWAP 


26.3 ±0.2 


29.6 ±0.2 


35.92 ±0.05 



Table 5: Cross-sections [fb] for e+e — > e z/ e r+iv 
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Figure 6: # e distribution and single-PF cross sections for ude v e . 
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Oe [Deg] 


FW 


EFL 


EFL/FW-1 (percent) 


0.0° 0.1° 


0.14154 


0.13448 


-4.99 


0.1° -^0.2° 


0.02113 


0.02031 


-3.88 


0.2° -^0.3° 


0.01238 


0.01194 


-3.55 


0.3° ^-0.4° 


0.00880 


0.00851 


-3.30 



Table 6: da/d9 e [pb/degrees], from WTO, for the process e + e — ► e v e u l _ i fi + , for |cos# e | > 0.997, 
Ep>15 GeV, and | cos 6» M | < 0.95. </s = 183 GcV. 





FW 


IFL 


IFL Q 


EFL 


EFL/FW-1 
(percent) 


183 GcV 


88.17(44) 


88.50(4) 


83.26(5) 


83.28(6) 


-5.5(5) 


189 GeV 


98.45(25) 


99.26(4) 


93.60(9) 


93.79(7) 


-4.7(3) 


200 GeV 


119.77(67) 


120.43(10) 


113.24(8) 


113.67(8) 


-5.1(5) 



Table 7: Cross section [fb] for the process e+e~ — > e~v e ud. M(ud) > 45 GeV, | cos0 e | > 0.997. No ISR. 





FW 


IFL 


IFL Q 


EFL 


EFL/FW-1 
(percent) 


183 GeV 


26.77(14) 


26.45(1) 


24.90(1) 


25.53(4) 


-4.6(5) 


189 GeV 


29.73(14) 


29.70(2) 


27.98(2) 


28.78(4) 


-3.2(5) 


200 GcV 


36.45(23) 


35.93(4) 


33.85(4) 


34.97(6) 


-4.1(6) 



Table 8: Cross section [fb] for the process e+e — ► e i/ e fj, + v /i . | cos# e | > 0.997, > 15GeV, | cos# M | < 
0.95. No ISR. 
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A second problem, relevant when including QED radiation, is the choice of the scales q 2 to be used 
in the Structure Function (SF) formalism, schematically represented in equation ([|) 

da = J dxi D{qj,Xi) da . (3) 

i 

The choice ~ s is proven to reproduce accurately the inclusive four-fermion cross sections, at least 
for s-channel dominated processes. For i-channel dominated processes, such as single-M^ production, a 
different choice is more adequate, as studies of certain processes have shown. When an exact first order 
QED radiative correction calculation exists for a t-channel dominated process, then the result can be 
compared to a Structure Function calculation with a q 2 scale related to the virtuality of the exchanged 
^-channel photon. With such a q 2 value the two kinds of calculations agree for small angle Bhabha 
scattering [[l9| and multi-peripheral two-photon processes [§0) , where the exact calculation already exists 
for some time |2^]. When no exact first order QED correction calculation is available, the first order 
soft correction may also serve as a guideline to determine q 2 , as proven by studies performed by GRACE 
and SWAP |2(], ^2|. In NEXTCALIBUR, the choice of the scale is performed automatically, event by event, 
according to the selected final state, as shown in table |)| The final state with 1 e~ (or 1 e + ) is relevant 
for single-W^ processes. 



Final State 


q 2 


2 


No e ± 


s 


s 


1 e - 


\t-\ 


s 


1 e+ 


s 


\t+\ 


1 e~ and 1 e + 


\t-\ 


\t+\ 


2 e~ and 2 e + 


min(|£_ |) 


min(|< + |) 



Table 9: The choice of the QED scale in NEXTCALIBUR. q± are the scales of the incoming while t± 
represent the i-channel invariants obtained by combining initial and final state e ± momenta. When two 
combinations are possible, as in the last entry of the table, that one with the minimum value of \t\ is 
chosen, event by event. 



A further problem is generating a p t spectrum for the photons. Usual solutions are the QED Parton 
Shower approach (QEDPS) and the use oip t dependent Structure Functions. 

For example pt dependent SF are implemented in NEXTCALIBUR via the replacement 



ml l-c + 2ff 9 2 (l-c + 2ff)* 



in the the soft part of the collinear SF, and 



(■?-) - I - 1 l-x 1 ( 



i_ c+2 ^ l + x 2 q 2 ( 1+a .)(x_ c + 2 ^)2 

2 

in the hard contributions. After integrating over c, one gets ln(- 2 -j) — 1. The inclusive QED result is 
therefore recovered and the photon spectrum is exact for small p t . 
Also GRACE and SWAP implement QEDPS and/or p t dependent SF. 

Since q 2 = s for s-channel legs and q 2 = \t\ for i-channcl dominated legs, the pt distribution is different 
in the two cases. In order to show this effect, the distribution in cos# 7 for the most energetic photon (with 
respect to the incoming e + ) in e + e~ — > e~v e ud (7), is shown in figure (?] as predicted by NEXTCALIBUR. 

The shape of the distribution can be easily understood in terms of radiation suppression. In fact 
\t\ — > implies -D(|i|, x) ~ 6(1 — x). Therefore the emission from the incoming e~ is suppressed, and the 
photons are preferably emitted collinearly by the incoming e + . 
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Figure 7: cos# 7 distribution for e + e 
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e~i> e ud(7), by NEXTCALIBUR. 



The question of comparing QEDPS and SF has been addressed by the GRACE group. The differences 
between the two approaches are at the order of 1% for single-14^ processes. Furthermore, the E~ distri- 
bution is well reproduced for soft photons, with respect to an exact Matrix Element calculations, while 
20% differences arc observed (and expected) in the hard region. 

In conclusion, putting together all the effects, a theoretical error of 5% is assigned to single- W 
processes. This result has to be compared with the accuracy required by the LEP2 collaborations namely 
2% for e + e~ — > e~ v e q q' (7), 5% for e + e~ — > eT v e e + v e (7), 5% for e + e~ — * e~ v e ft + (7) 0. 
However, different pieces of knowledge are still scattered in different codes, and an improvement of the 
present situation is possible via a multi-step experimental procedure, in which, for example, code A is 
used to correct code B for the missing effects. 

At the LC further improvements are needed. In particular, the genuine EW corrections get larger with 
increasing energy. Again, a full one-loop four-fermion calculation seems unavoidable, for high precision 
measurements . 

3 Four fermions plus 1 visible photon 

This signature gives information on the quartic gauge coupling (see figure and is relevant when studying 
processes with three final state bosons, such as W + W~j production, ZZj and Z77. Furthermore, it is 
a building block for the full computation of e + e~ — > 4f at 0(a). 

A bunch of codes contributed, with different strategies. CompHEP, GRACE and HELAC/PHEGAS compute 
the exact Matrix Element (ME) with massive fermions. RACDDNWW uses exact ME, but in the limit of 
massless fermions. NEXTCALIBUR generates photons only through pt dependent SF. WRAP has a matching 
between ME and SF generated photons. The last approach allows to estimate the size of the double 
counting when blindly dressing the 4/ + 7 ME with collinear ISR. WRAP observed effects up to 5%, 
depending on the energy cut used to define the visible photon. 

2 For semi-leptonic processes this number refers to the High Mass region, namely to situations when appropriate cuts on 
the jet-jet system are applied in order to suppress the non-perturbative resolved photon contributions. 



In figure |9j a study by CompHEP is shown on the reliability of a narrow-width approximation. 
Distributions generated by GRACE are shown in figure [To| . 




Figure 8: Quartic gauge boson coupling in four-fermion production plus 1 visible 7. 




Figure 9: Distributions in the 7 energy, 7 transverse momentum, 7 angle with the beam, and in the 
opening angle between the 7 and the nearest charged fermion. The distributions for the e + e _ — > ^fiP^ud 
are shown by the solid line and the distributions for the e + e" — > r yfiP Pl W + with the following W isotropic 
decay are shown by the dashed line. 
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Figure 10: E 11 COS Q~y and invariant mass distributions for [iv^ud^. 
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In tables |o| and |Tl| the cross sections for inclusive (crt t)j non radiative (cr nra d) , single radiative (<7 sra d) 
and double radiative (crdrad) events are shown, as computed by NEXTCALIBUR. This kind of study is useful 
for signal definition. 



Type 


Cross-section 


Ctot 


16.107(9) 


^nrad 


15.018(9) 


Csrad 


1.0697(30) 


fdrad 


0.0189(4) 



Table 10: a in fb for e + (l)e"(2) \T (3) fi+ (4)u(5)u(6) . M(34) > lOGeV and M(56) > lOGeV. ZZ 
like cuts, E y > 1 GeV, | cos6» 7 | < 0.985. 



Type 


Cross-section 


Otot 


617.27(59) 


^nrad 


578.19(58) 


^srad 


38.54(16) 


Cdrad 


0.54(2) 



Table 11: a in fb for e+(l)e-(2) -> /*-(3)z^(4)«(5)d(6). M(56) > lOGeV. V^W^ like cuts, £; 7 > 
lGeV, |cos6» 7 | < 0.985. 

Finally, in table [l2| and figure O, I show results of tuned comparisons among HELAC/PHEGAS, RAC00NWW 
and WRAP. 

In conclusion, a very good technical precision has been reached in the computation of four-fermion 
processes plus 1 additional photon. However, the non-logarithmic 0(a) corrections are not known. 
Therefore a 2.5% theoretical accuracy on total cross section and inclusive distributions is estimated at 
LEP2 energies |Q. Larger effects are expected at the LC. 

4 Z-pair production 

Less accuracy is required at LEP2 for this observable with respect to the VF-pair case. Thought feasible in 
principle, a DPA ZZ calculation is not available yet. A theoretical accuracy of 2% on azz is estimated at 
LEP2 by varying the renormalization scheme and by comparing different codes and different treatments 
of the QED radiation (see table [l3]) . 
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UR AP 




PHFHA^ /HFT Af 
± Oi_ivji-i.o / niLJ_ifio 


ude~v e "f 
udscq 


75.732(22) 
78.249(43) 
28.263(9) 
29.304(19) 
199.63(10) 


75.647(44) 
78.224(47) 
28.266(17) 
29.276(17) 
199.60(11) 


75.683(66) 
78.186(76) 
28.296(22) 
29.309(25) 
199.75(16) 



Table 12: a in fb from WRAP, RacoonWW and PHEGAS/HELAC. 



1.03 




-1 -0.5 0.5 1 
COS6 y 

Figure 11: cos# 7 distributions and ratios for v^^ + e"v e ^i. 



channel 


YFSZZ 


ZZTO GF-scheme 


ZZTO a-schcmc 


qqqq 


294.6794(490) 


298.4411(60) 


294.5715(59) 


qqw 


175.4404(302) 


175.5622(35) 


174.9855(35) 


qq\\ 


88.1805(134) 


88.7146(18) 


87.9881(18) 


\\w 


26.2530(463) 


26.0940(5) 


26.1342(5) 


1111 


6.5983(15) 


6.5929(1) 


6.5706(1) 


WW 


26.1080(71) 


25.8192(5) 


25.9868(5) 


total 


617.2596(755) 


621.2241(124) 


616.2366(123) 



Table 13: ZZ cross-section [fb] at y/s = 188.6 GeV. YFS exponentiation is used by YFSZZ, while ZZTO 
includes EFL corrections in the a-scheme. 
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5 Conclusions 



Four-fermion Physics at LEP2 is in a good shape. Thanks to the results of ref. |l|, improved calculations 
are available for 

• W-pair production 

• Single-VF production 

• Four fermions plus 1 visible 7 

• Z-p&ir production. 

In general, the theoretical accuracies required by the LEP2 experiments are achieved. 

At the LC more accuracy is needed. In particular radiative corrections must be included. 

Four-fermion loop calculations at 0{a s ) already exist plj] , while a full 0(a) EW calculation is still 
missing. 
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